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ABSTRACT 

An  extension  of  two-color  Particle  Image  Velocime- 
try  (PrV)  is  described  in  which  the  recording  media  (color 
film)  is  replaced  with  a  high-resolution  (3060  x  2036  pixel) 
color  CCD  sensor.  Incorporation  of  the  sensor  eliminates 
the  time  associated  with  film  development  and  digitization. 
Use  of  color  allows  the  removal  of  directional  ambiguities 
without  resorting  to  polarization-based  image-shifting  tech¬ 
niques.  For  comparing  the  performance  of  the  color  CCD 
sensor  with  conventional  color  film,  PIV  images  were 
obtained  on  a  sand-blasted  surface  undergoing  a  known 
translation  and  rotation  under  comparable  lighting  and  cam¬ 
era  conditions.  Good  agreement  was  obtained  between  the 
two  recording  media,  indicating  that  the  behavior  of  the 
color  CCD  sensor  is  similar  to  that  of  400-ASA  color  film. 
This  digital  two-color  PIV  (D2CPIV)  technique  was  used  for 
studying  simulated  turbine  film-cooling  flows  to  obtain 
more  detailed  characterization  of  the  coolant-injection  phe¬ 
nomena  and  their  interaction  with  freestream  disturbances. 
This  technique  allowed  near-real-time  optimization  of  the 
PIV  parameters  which  resulted  in  higher  valid  vector  density 
and  shear-layer  resolution. 


NOMENCLATURE 

d  film-cooling-hole  diameter  (1.905  cm) 

R  coolant  blowing  (mass  flux)  ratio  (p,U/p^UJ 
Re  Reynolds  number  based  on  film-cooling-hole 
diameter 

Tu  turbulence  intensity  (u/U) 

U  mean  local  streamwise  velocity  (m/s) 

X  streamwise  distance  measured  from  the  downstream 
lip  of  the  injection  hole  (cm) 
y  vertical  distance  from  the  injection  surface  (cm) 


1.  INTRODUCTION 

The  PIV  technique  has  been  in  use  for  a  number  of 
years  to  measure  velocity  distributions  in  planar  cross  sec¬ 
tions  of  aerodynamic  flowfields  (Adrian  1991).  One  of  the 
difficulties  involved  in  implementing  this  velocimetry 
technique  is  the  180-deg.  directional  ambiguity  which 
results  from  the  inability  to  determine  the  temporal  sequence 
of  the  particle  pairs.  Several  techniques  have  been  devel¬ 
oped  to  resolve  this  ambiguity  problem;  most  involve 
imposing  a  shift  between  consecutive  image  exposures  by 
means  of  scanning  or  rotating  mirrors  (Adrian  1986),  pulse 
tagging  (Grant  and  Liu  1990),  calcite  crystals  (Landreth  and 
Adrian  1988),  or  polarizing  beam  splitters  (Lourenco  1993). 
For  overcoming  the  difficulties  inherent  in  these  techniques, 
a  two-color  PIV  system  was  developed  (Goss  et  al.  1991). 
The  advantages  of  this  system  are:  1 )  the  directional  ambi¬ 
guity  is  resolved  using  the  color  coding  of  the  particle 
images,  which  is  inherent  in  the  system,  2)  higher  data 
yields  and  signal-to-noise  levels  are  attainable,  and  3)  the 
technique  is  suitable  for  both  reacting  and  non-reacting 
flowfields. 

In  early  experimental  approaches  involving  the  PIV 
technique,  the  particle  images  were  recorded  on  photo¬ 
graphic  film.  However,  this  type  of  recording  is  time  con¬ 
suming  because  of  the  need  to  develop  the  film  before 
digitization  and  subsequent  computer  processing.  This  dis¬ 
advantage  can  be  overcome  by  recording  the  particle  images 
directly  onto  a  two-dimensional  CCD  array.  This  approach 
has  been  recognized  by  several  investigators.  Cho  0989) 
proposed  that  the  double-exposed  digital  images  be  obtained 
by  digitizing  single-exposed  video  images  and  adding  the 
successive  images.  Okada  et  al.  (1990)  used  liquid-crystal 
television  (LCTV)  having  a  resolution  of  320  x  320  pixels. 
Willert  and  Gharib  (1991)  developed  digital  PFV  and  used 
cross  correlation  of  the  original- image  frames  acquired 
in  succession  to  obtain  an  accurate  measurement  of  the 
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displacements  without  directional  ambiguity.  Lourenco  et 
al.  (1994)  introduced  a  fully  digital  and  operator-interactive 
PIV  system  which  utilized  a  CCD  area  sensor  (1320  x  1035 
pixel),  resulting  in  significant  improvements  over  earlier 
CCD-based  PIV  systems.  This  system  has  been  successfully 
implemented  in  a  wide  variety  of  flow  regimes. 

The  extension  of  the  two-color  PIV  technique  to 
include  CCD  cameras  has  been  hampered  in  the  past  by  the 
lack  of  commercially  available  high-resolution  color  CCD 
camera  systems.  Goss  et  al.  (1994)  attempted  to  utilize  a 
single  CCD  camera  and  dichroic  mirror  to  record  the  two- 
color  PIV  images  of  a  stagnation-point  flow.  However, 
several  problems  were  encountered,  including  a  nonuniform 
shift  across  the  image  and  degraded  signal-to-noise  ratios. 
Because  of  recent  developments  in  high-resolution  color 
cameras,  the  difficulty  in  utilizing  color  CCD  cameras  for 
two-color  PIV  has  significantly  decreased. 

The  present  paper  describes  the  extension  of  two- 
color  PIV  by  recording  the  color  images  onto  a  single,  high- 
resolution,  digital  (3060  x  2036  pixel)  color  CCD  sensor, 
thus  eliminating  the  processing  time  and  subsequent  digiti¬ 
zation  time  of  color  film  and  the  complexities  associated 
with  conventional  image-shifting  techniques.  For  demon¬ 
strating  the  direction-resolving  capabilities  of  the  D2CPIV 
technique  and  evaluating  the  spatial  resolution  capabilities, 
a  test  was  performed  using  a  sand-blasted  rotating  wheel. 
Results  from  this  study  indicate  that  the  high-resolution 
color  CCD  sensor  is  comparable  to  400-ASA  color  film  in 
sensitivity,  spatial  resolution,  and  data  quality. 

This  D2CPIV  technique  was  also  used  to  study  simu¬ 
lated  turbine  film-cooling  flows  in  order  to  provide  a  more 
detailed  characterization  of  the  coolant-injection  phenom¬ 
ena  and  their  interaction  with  the  freestream  disturbances. 
These  types  of  flows  are  important  because  many  state-of- 
the-art  turbine  stages  employ  film  cooling  to  permit  near- 
stoichiometric  combustor  operating  temperatures.  Film¬ 
cooling  air  is  injected  through  rows  of  small  (0.5  -  0.8  mm 
diameter  typical)  holes  in  the  blade  surface.  The  coolant  air 
is  supplied  from  the  compressor  exit  flow  and  is  maintained 
at  essentially  constant  pressure. 


2.  EXPERIMENTAL  SETUP  AND  PROCEDURES 
2.1  Facility  Description 

The  open-loop  film-cooling  wind  tunnel,  as  shown 
in  Figure  1,  has  been  described  in  detail  by  Bons  et  al. 
(1994,  1995).  The  main  flow  passes  through  a  conditioning 
plenum  containing  perforated  plates,  honeycomb,  screens, 
and  a  circular-to-rectangular  transition  nozzle.  Downstream 
of  the  transition  nozzle,  at  the  location  of  the  film-cooling 
station,  low  freestream  turbulence  levels  of  0.7%  (±  0.05) 
can  be  achieved,  with  velocity  and  temperature  profiles 
being  uniform  to  within  1%  at  the  film-cooling  station. 
High  freestream  turbulence  levels  up  to  17%  can  be  achieved 
at  the  film-cooling  station.  A  single  row  of  1.905-cm  film¬ 
cooling  holes  at  an  injection  angle  of  35  deg  to  the  primary 
flow  is  investigated.  The  length-to-diameter  ratio  of  the 


film-cooling  holes  evaluated  is  2.4.  The  ratio  of  integral 
turbulence  scale  to  film-hole  diameter  is  in  the  range  2.88  - 
4.23,  depending  on  the  turbulence  level  and  turbulence- 
generation  mechanism  (Bons  et  al.  1994),  The  ratio  of 
momentum  thickness  to  hole  diameter  typically  is  0.05. 
The  film-cooling-hole  Reynolds  numbers  (Re)  are  10,000, 
20,000,  and  40,000.  The  ratio  of  micro  scale  to  film-hole 
diameter  is  in  the  range  0.1  -  0.39.  The  ratio  of  temperature 
or  density  of  the  film  flow  to  the  primary  flow  is  typically 
in  the  range  1.07  -  1.09.  The  variation  in  blowing  ratios  (R 
=  0.7  -  1.5)  is  achieved  effectively  by  varying  the  velocity 
ratio  of  the  film-flow  to  the  primary-flow  velocity. 

2.2  Two-Color  PIV  System 

The  two-color  PIV  system  uses  color  for  temporal 
marking  of  the  seed  particles  in  the  flow  field.  The  green 
(532-nm)  laser  output  from  a  frequency-doubled  Nd:YAG 
laser  and  the  red  (640-nm)  laser  output  from  a  Nd:YAG- 
pumped  dye  laser  (DCM  dye)  are  combined  by  a  dichroic 
beam  splitter  and  directed  through  sheet-forming  optics. 
The  laser-sheet  energy  is  typically  20  mJ/pulse,  with  a 
thickness  of  <  1  mm  at  the  test  section.  The  temporal  delay 
between  the  two  lasers  is  a  function  of  gas  velocity,  optical 
magnification,  and  interrogation  spot  size.  In  the  present 


Fig.  I  Schematic  diagram  of  experimental  setup,  (a)  side 
view,  (b)  top  view. 
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experiments,  the  time  delay  between  the  two  color  lasers  is 
set  at  10  psec  for  flows  at  Re  =  40,000. 

The  film-cooling  flow  is  seeded  with  sub-micron  size 
particles,  and  the  Mie  scattering  from  the  seed  particles  is 
recorded  on  a  Kodak  DCS  460  CCD  array.  This  CCD  sensor 
has  a  resolution  of  3060  x  2036  pixels,  and  each  pixel  is  9 
pm  square.  The  CCD  camera  has  a  built-in  12-bit  analog-to- 
digital  converter  for  increased  dynamic  range  and  a  frame 
rate  of  I  frame/sec.  It  also  features  a  PC-MCIA  storage  drive 
which  delivers  about  26  exposures,  with  each  PIV  image 
occupying  -  18  Mbytes.  A  105-mm  micro  lens  with  an  f- 
number  5.6  is  used  to  record  the  images. 

The  color  response  of  the  DCS  460  CCD  array  is 
accomplished  by  overlaying  the  individual  pixels  of  the 
camera  with  a  series  of  red,  green,  and  blue  color  filters.  The 
green-  and  red-laser  outputs  are  situated  near  the  peak  of  the 
transmission  of  the  green  and  red  camera  filters,  respec¬ 
tively.  Because  the  camera  was  built  to  respond,  in  part,  as  a 
human  eye  to  color,  most  of  the  pixels  are  green-sensitive. 
The  relative  percentages  of  green,  blue,  and  red  pixels  are 
50%,  25%,  and  25%  respectively.  The  distribution  of  the  red 
and  blue  pixels  is  random  for  minimizing  straight-edge 
effects.  The  output  from  the  camera  controller  to  the  com¬ 
puter  is  a  24-bit  RGB  tiff  image;  thus,  the  12-bit  ADC 
camera  output  must  be  mapped  into  three  8-bit  colors,  each 
having  a  spatial  resolution  equivalent  to  the  chip  size  (3060 
X  2036).  This  is  accomplished  through  proprietary  software 
developed  by  Kodak  which  involves  interpolation  to 
increase  the  spatial  resolution  of  the  camera.  Because  of  the 
proprietary  nature  of  this  interpolation  software,  it  was  not 
known  whether  the  resulting  resolution  of  the  color  camera 
would  be  sufficient  for  accurate  PIV  measurements.  This 
study  was  also  directed  toward  evaluation  and  demonstration 
of  the  use  of  the  Kodak  color  CCD  camera  for  PIV  studies. 

2.3  Data  Analysis 

Once  the  PIV  image  has  been  capmred  and  digitized, 
the  velocity  field  is  obtained  using  a  cross-correlation  tech¬ 
nique.  To  improve  the  analysis  of  the  seeded  flow  field,  the 
output  of  the  linear  camera  was  convolved  with  a  logarithm¬ 
like  function  prior  to  cross-correlation  analysis.  The  pres¬ 
ent  cross-correlation  technique  is  based  on  intensity  maps 
of  the  red  and  green  images  of  the  scattered  light. 

Consider  the  intensity  distributions  of  the  red  and 
green  images  r{x,y)  and  g(x,y)  and  their  corresponding 
Fourier  transforms  R(a,p)  and  G(a,p).  The  two-dimensional 
cross-correlation  function 

h(x,y)  =  Jr  Jr  r(a.p)  g(x+a,y+P)  da  dp 
=  F  [F(r(x.y))  F(g'(x,y))] 

=  f' [R(a.p)  G’(a,P)]  (D 

is  employed  to  determine  the  magnitude  and  direction  of  the 
average  velocity  over  the  interrogation  area.  (Note  that 
unlike  in  processing  methodologies  that  are  based  on 
autocorrelation,  the  direction  of  the  velocity  vectors  is 
uniquely  determined.) 


The  correlation  function  is  calculated  over  small 
segments  (interrogation  domains)  of  the  PIV  image.  The 
dimensions  of  each  interrogation  domain  are  dependent  on 
particle  density,  estimated  local  velocity  gradients,  particle- 
image  size,  and  desired  spatial  resolution.  The  maximum 
displacement  of  each  particle  must  be  less  than  one-half  the 
interrogation  spot.  In  the  present  experiments,  the  interro¬ 
gation  domain  measured  64  x  64  pixels,  corresponding  to  2 
X  2  mm  in  the  measured  flow.  For  enhancing  the  overall 
resolution,  the  interrogation  domains  are  overlapped  by 
one-half  the  domain  size.  The  peak  of  the  correlation  map 
corresponds  to  the  average  velocity  displacement  within  the 
interrogation  spot.  An  intensity- weighted  peak-searching 
routine  is  used  to  determine  the  exact  location  of  the  peak  to 
sub-pixel  accuracy.  The  number  of  particle  pairs  normally 
necessary  to  ensure  a  desirable  signal-to-noise  ratio  is 
reduced  to  four  or  five  pairs  when  the  cross-correlation 
analysis  is  employed. 

2.4  Uncertainty  Analysis 

The  experimental  uncertainties  are  calculated  based 
on  knowledge  of  the  instrumentation  used  and  a  simple  root- 
mean-square  error  analysis  (Kline  and  McClintock  1953). 
This  method  assumes  that  contributions  to  uncertainties 
arise  mainly  from  unbiased  and  random  sources.  Uncertainty 
in  the  velocity  measurement  arises  from  the  time  required  to 
keep  the  large  out-of-plane  velocities  and  fluctuating  com¬ 
ponents  within  the  laser  sheet  during  both  pulses.  The 
resulting  number  of  pixels,  typically  10,  of  displacement 
and  the  sub-pixel  resolution  of  0. 1  pixels  then  dictate  the 
uncertainty  of  ±  1%.  The  spatial  resolution  in  the  data  pre¬ 
sented  is  4  mm.  The  data  are  acquired  at  a  resolution  of  2.0 
mm  using  32  pixels/mm  for  a  spatial-resolution  accuracy  of 
-  0.3%. 


3.  RESULTS  AND  DISCUSSION 
3.1  Calibration  of  the  System 

For  testing  the  ability  of  the  color  CCD  camera  to 
capture  images  of  sufficient  quality  for  PIV  analysis  and  for 
demonstrating  the  directional  resolving  ability  of  the 
two-color  PIV  system,  “calibration”  displacements  were 
obtained  by  recording  the  speckle  patterns  from  a  sand¬ 
blasted  surface  undergoing  uniform  translation  and  rotation. 
The  output  from  the  green  and  red  lasers  was  sent  unexpanded 
onto  a  scattering  plate  which  diffusely  scattered  the  light 
onto  the  sand-blasted  surface.  In  the  case  of  uniform  transla¬ 
tion,  the  displacements  for  each  translation  were  generated 
by  recording  the  single-exposed  digital  image  (green  at  the 
initial  position  and  red  after  a  small  translation)  and  then 
numerically  adding  the  two  images.  Extreme  care  was  taken 
to  prevent  any  movement  of  the  camera  or  the  stage  between 
exposures.  Cross-correlation  analysis  was  employed  to 
process  these  images  using  a  64  x  64  pixel  interrogation 
size.  The  processed  translation  data  which  contain  about 
6000  vectors  showed  uniform  displacements  throughout  the 
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interrogation  region,  and  the  rms  value  of  the  translational 
horizontal  component  was  found  to  be  0.01 1  mm.  The  dis¬ 
placements  for  the  uniform  rotational  case  were  generated  by 
spinning  the  wheel  and  exposing  the  speckle  pattern  from 
both  green  and  red  lasers  on  the  same  frame.  Figure  2(a) 
shows  a  typical  speckle  pattern  of  the  rotational  image,  and 
Figure  2(b)  shows  the  velocity  data  which  clearly  indicate 
the  ability  of  the  two-color  approach  to  resolve  the  direc¬ 
tional  ambiguity.  In  order  to  compare  results  obtained  with 
the  color  CCD  sensor  and  film,  a  line  passing  through  the 
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Fig.  2  Rotational  data  showing  the  direction-resolving 
capability  of  the  D2CPIV  technique,  (a)  speckle 
pattern,  (b)  velocity  field 


center  of  rotation  with  zero  horizontal  displacement  was 
examined.  Because  the  tangential  velocity  of  a  rigid  body 
undergoing  rotation  varies  directly  with  the  radial  distance 
from  its  center  and  its  rotational  speed,  a  simple  deter¬ 
ministic  relationship  exists  between  the  actual  and  measured 
velocities.  Velocity  measurements  were  obtained  using  both 
the  color  CCD  sensor  and  color  film  (400  ASA)  under  identi¬ 
cal  conditions  such  as  magnification,  recording  lens,  aper¬ 
ture,  and  rotational  speed.  The  film  image  was  digitized  at  a 
resolution  of  2700  pixels/in.,  yielding  an  image  size  of 
-  3000  X  2000  pixels. 

These  images  were  interrogated  using  spot  sizes 
ranging  from  32  x  32  to  256  x  256  pixels.  Figure  3  shows  a 
comparison  of  the  actual  and  measured  velocities  for  both 
film  and  CCD  sensor,  with  the  images  being  processed  using 
a  128  X  128  spot  size.  This  shows  that  the  data  obtained 
using  the  high-resolution  CCD  sensor  are  very  comparable 
to  the  film  data.  The  proprietary  interpolation  software 
developed  by  Kodak  appears  to  be  effective  in  maintaining 
equal  red,  green,  and  blue  pixel  spatial  resolution,  even 
though  the  color  pixels  are  not  distributed  equally  through 
the  color  CCD  sensor. 

3.2  Application  to  Simulated  Turbine 
Film-Cooling  Flows 

For  further  demonstrating  the  ability  of  the  system  to 
capture  the  complex  flow  features  in  quantitative  form, 
simulated  turbine  film-cooling  flows  were  considered.  These 
flows  occur  in  a  very  hostile,  unsteady  environment  where 
velocity  and  temperature  disturbances  exceed  20%.  The 
freestream  and  film-cooling  wall  conditions  in  these  flows 
exceed  those  encountered  in  a  classical,  fully  turbulent 
boundary  layer.  Accurate  modeling  of  these  flows  has 
proved  to  be  difficult  due  to  the  high-level  unsteadiness. 
These  flows  have  been  investigated  by  Gogineni  et  al. 


Fig.  3  Comparison  of  film  and  CCD  data 
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(1995,  1996)  using  the  conventional  two-color  PIV  tech¬ 
nique  for  Re  =  20,000,  with  R  varying  from  0.5  to  1.5  and 
Tu  from  1  to  17%.  It  was  shown  that  turbulence  clearly 
increases  the  film  spread  for  R  =  0.7.  It  was  also  observed 
that  a  linear  relationship  exists  between  the  jet-exit  slopes 
and  both  Tu  and  R  for  R  <  1 .0. 

In  the  present  study  PIV  measurements  were  made 
using  a  CCD  sensor  for  Re  =  10,000,  20,000,  and  40,000  at 
freestream  Tu  levels  of  1  and  17%  by  varying  R  from  0.7  to 
1.5.  These  flow  conditions  were  selected  to  understand  the 
influence  of  high  Reynolds  number  and  high  freestream 
turbulence  on  film-cooling  effectiveness.  Since  the  instan¬ 
taneous  realizations  are  not  representative  of  mean  flow 
behavior,  10  images  were  recorded  for  each  condition.  Fig¬ 
ures  4(a)  -  4(c)  show  typical  double-exposed,  two-color  PIV 
images  for  R  =  0.7  and  Tu  =  1%  as  Re  is  varied  from  10,000 
to  40,000.  In  these  photographs,  only  the  film-cooling 
flow  is  seeded  with  sub-micron  size  particles.  When  Re  is 
increased,  the  shear-layer  frequency  over  the  film  hole 
increases  dramatically.  The  dominance  of  the  film  hole  in 
setting  the  shear  layer  frequency  over  the  film  hole  is 
reduced  to  0.75  of  the  film  hole  diameter  as  Re  is  increased 
to  40,000.  In  ail  three  photographs,  a  second,  parallel 
shear  layer  is  detected  in  the  leading-edge  exit  region.  This 
can  be  associated  with  separation  from  the  inside  entrance 
lip  of  the  film-cooling  hole.  The  shear-layer  roll  up  is 
observed  to  be  opposite  that  which  would  be  expected  for 
R  =  0.7.  The  observed  roll  up  results  from  the  film  tube 
boundary  layer  and  a  lack  of  seeding  in  the  freestream.  Fig¬ 
ures  4(d)  -  4(f)  show  the  corresponding  instantaneous 
velocity  distributions.  For  clarity  purposes,  only  a  few  of 
the  measured  vectors  are  shown  in  each  frame. 

Figures  5(a)  -  5(c)  show  the  PIV  images  when  R  is 
increased  to  1.5  and  Tu  to  17%.  At  R  :=  1.5,  the  scale  of  the 
shear  layer  after  the  film-cooIing  hole  increases  or  the  fre¬ 
quency  decreases  as  compared  to  R  =  0.7  in  Figure  4.  Figures 
5(d)  -  5(f)  show  the  corresponding  instantaneous  velocity 
dis-tributions.  The  CCD  digitization  allowed  near-real-time 
optimization  of  the  seeding  and  laser-sheet  intensities 
which  resulted  in  higher  valid  vector  density  and  shear-layer 
resolution. 


4.  CONCLUSIONS 

A  new  implementation  of  the  two-color  PFV  tech¬ 
nique  is  described  where  conventional  color  film  is  replaced 
with  a  high-resolution  color  CCD  sensor.  The  proprietary 
interpolation  software  developed  by  Kodak  is  effective  in 
maintaining  equal  red,  green,  and  blue  pixel  spatial  reso¬ 
lution,  even  though  the  pixels  are  not  equally  distributed.  It 
was  determined  that  optimum  results  in  seeded  flow  fields  are 
obtained  by  applying  a  logarithm-like  function  to  the  cam¬ 
era  image  prior  to  cross-correlation  analysis.  Tests  on  a 
rotating  wheel  showed  the  directional-resolving  capabilities 
of  the  two-color  approach,  and  the  system  was  successfully 
applied  under  simulated  turbine-flow  conditions.  The  CCD 


sensor  contributed  significant  improvement  to  the  optimi¬ 
zation  of  seeding  and  laser-sheet  intensities  for  the  simu¬ 
lated  turbine  high  freestream  flows. 
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